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Abstract 

We calculate the nuclear suppression for J/ty and VP' production within a coupled 
channel approach in the subspace of the J/vl/ and VP' states. We are able to explain, 
why (i) the J /vP and VP' show the same suppression from 200 GeV to 800 GeV in 
proton- nucleus collisions and why (ii) the is absorbed more strongly than the J/fy 
in nucleus-nucleus collisions at 200 GeV. The numerical result which includes only 
interactions with nucleons acconts for half of the observed suppression in sulphur- 
uranium collisions. 



The E772 collaboration |L| was the first who claimed that $ (we use hereafter this 
abbreviation for J/^) and ^' produced in proton-nucleus (p — A) collisions at 800 GeV and 
x F > 0.1 experience the same nuclear suppression. This result was confirmed by the NA38 
collaboration 0, |3|, |4j] at an energy of 200 GeV and for ij? « 0.1, though with larger error 
bars. The two observations contradict the simple-minded expectation, namely that the VP' 
should be more strongly suppressed, because absorption cross section scales with the mean 
square radius (r 2 ) of the meson |5|, || and the ^' is much larger than the To our opinion, 
no satisfactory explanation of the data has yet been proposed. 

Recently the situation became even more mysterious by the observation & Q that in 
Sulfur-Uranium (S-U) collisions at 200 GeV /A the ^' is significantly more strongly absorbed 
than the Does this result signal dense hadronic gas or quark-gluon plasma formation? 

In this letter we treat the propagation of a cc pair through nuclear matter as a coupled 
system of the \1> and ^' states. In addition to elastic collisions \I>iVfi\I>iV and ^'NM'N 
with amplitudes f(^, and f(^', \&') respectively, we consider the conversion amplitudes 
f(^, and f(^f', ^) for the processes # ty'N during propagation through nuclear 
matter. The inelastic amplitudes turn out to be nearly as big as the elastic ones. 

We define as nuclear suppression factor of a \1/ meson in pA collisions the ratio 

A ^ _ 1 a(pAmX;x F} E)/dx F 

* [XF: ] ~ A a( P NWX; x F , E)/dx F ' 1 j 

where E is the lab. energy. We also introduce the relative to \I/ nuclear suppression 
function by 

qpA _ qpA I qpA I cy\ 



Then the results of the E772 and NA38 experiments can be summarized by S^f/^ ~ 1 



for the values of x F and E considered. 

In order to exhibit the physics of our coupled channel approach, we first calculate eq. 
(Q) perturbatively, i.e. restricting ourselves to only one fif'N) interaction after the 
creation of the cc pair. This interaction can be an elastic one as well as a conversion event. 
Then 



qP A ^-\°t?(r + e/R){T) A 
V/9 l-\af ?{l + eR){T) A ' 1 j 

where r = /(\^'\I>')//(\&^), e = / '{^') f and cr^f = 2Imf(^, \I>). Furthermore, 
i? denotes the relative amplitude of ^ to \I/ production in the initial pN collision, and 
(T)a = 1/A J d 2 bT 2 (b), where T{b) = J^dzpA^b, z) is the nuclear thickness. 

If one neglects the conversion rate (e = 0), eq. (^) reduces to the conventional result, 
namely that the \I/ suppression depends only on af a f , while raf a f = af^ t N determines the 
suppression of the 

The ratios r and e in eq. (0) can be fairly reliably calculated since at high energies 
the scattering amplitudes f(a,/3), where a and (3 stand for \& and are proportional 
to (a\r^\/3), as long as the meson radius in the transverse direction, (rf,), is small. With 
IS* and 25' harmonic oscillation functions for \1/ and respectively, one has r = 7/3 and 



e = — y 2/3. The ratio R of initial \l/ to \&' production cannot be calculated in this way. We 
turn the argument around and calculate from eq. ([![) that value of R ca i which leads to the 
observed relative suppression S^,,^ = 1 and compare it to R exp deduced from pN collisions. 



A quadratic equation for R leads to R ca i = y5/3±yz/3, the smaller one being R ca i = 0.47. 
If one uses the experimental intensities of and \I/ produced in pp collisions and corrects 
them for the feedings ^I/'fi^ and i one arrives at an amplitude ratio \R exp \ = 0.48 ±0.06 
]7|. f| , which agrees well with the calculated one, indicating that the initially produced state 



\&i C ) = (l^) + R\^'))/y \ + R? is such as to lead to the same final state attenuation for \l/ 
and ^' . Mathematically spoken, |$^ c ) is an eigenstate of the final state interaction matrix, 
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/(*,*)• (4) 



The property of |$f c ) being eigenstate of / is equivalent to the statement that |$^ c ) has 
an extreme value for its transverse size: 

A ( $ r -| r 2|$f ) = . (5) 
The experimental value \R eX p\ selects the physical state as that with minimal transverse 



extension. 

Strictly speaking the result eq. (|3|) is only valid for energies EQoo. For finite lab. 
energies, especially for pA and AA collisions at 200 GeV one has to include the effect of the 
longitudinal momentum transfer q associated with the conversion reaction \l/ # ty': 

Ml - Ml 

q= -2E^T> (6) 



where x\ = (xf + yxp + 4M|/s)/2 with yfs the cm. energy. One arrives at an expression 
like eq. @ where e is replaced by eF A (q) with 

2 poo poo 

F A(q) = -jj^r J_^dzp A (b,z) dz' p A (b, z') exp(iqz') (7) 

being a kind of nuclear formfactor. For the E772 experiment q < 0.06 frrT 1 for xp > 0.2 and 
Fa = 1 is a very good approximation. Thus, the observed result S^,,^ = 1 is reproduced. 

For the NA38 experiment at 200 GeV/ A especially for the nucleus-nucleus collisions, 
the effect of the formfactor becomes crucial: In the S — U collision the produced $ (or 
moves with fractional momentum xp (in the cm. system) with respect to the target nucleus 
U, but moves with — xp relative to the projectile nucleus S (inverse kinematics). Therefore 
the suppression arises from two sources, 

7* = »S#,^(— xp, E)Sq,,/y Sf (xF, E) . (8) 

For 200 GeV and xp = 0.2 we have q = 0.16 jm~ x to be used in the second factor of eq. (^) 
and q = 0.56 /m _1 in the first one. This is the inverse kinematics which leads to a much 
stronger ty'/ty suppression. 

For the detailed comparison with experiment we have to go beyond the perturbative 
expression @ and use numerical methods. We describe the propagation of cc pair through 
nuclear matter by a differential equation || |J 

i^ c %z,b)) = U(z,b)\^(z,b)) , (9) 

where the limitation to the ^, ^' subspace invites the use of matrix notations: |$ cc ) = f^J 
and 



i \ ° ° -k^V(M) 1 e ! do) 

q J A \ e r 

with initial condition |$^) = ( jj) / Vl + R 2 at the point (6, zq) of cc creation. At z — +oo 
the wave function |$ cc ) is projected on the \l/ and states. The result is squared and 
averaged over the coordinates (b, z ) of the production point. Note that this is only true if 
the longitudinal momentum transfer at the production point, go = M\j2Ex\ ^> q, is much 
larger than the reversed mean internucleon distance in a nucleus. Otherwise one should take 
into account coherence between different production points as well (see discussion in terms 



of production and formation times in [HJ 11]). This would be equivalent to an effective 
increase of the length of path of the qq pair in nuclear matter. However it does not affect 
the relative production rate if |$ cc ) is an eigenstate of interaction. 

The relative suppressions are calculated for p — A collisions and with the help of eq. 
(H) also for A — A collisions. The numerical results shown in Figs. 1 and 2 are calculated with 
realistic nuclear densities and with the values of e and r as given by the harmonic oscillator 
model and for the absolute value of the total cross section afj* ~ C^rf,)^ = 5.7 mb, 



where we use the perturbative QCD estimate of [1 1] or systematics of || for a value of C. 
While at 800 GeV the values for S& /9 are measured directly, the values given at 200 GeV 
for B^/a^i / B^a^ in nuclear collisions were renormalized by us using the value 1.80 ±0.10 % 
for this quantity from pp and pd collisions 0] . Fig. 1 shows the suppression functions SV m 
for p — W at 800 GeV and for p — U and p — W at 200 GeV. Although we predict for 
200 GeV some reduction of the ty'/ty relative suppression it is still in agreement with the 
data of the NA38 experiment within rather large error bars. Fig. 2 shows the predicted and 
observed |fj, |] suppression 5*f for the nucleus-nucleus case. The solid curve is the product 
of the suppression curves for the p — U and the S — p collisions (each dashed). Experiment 
and calculation agree in that the should be significantly more strongly suppressed than 
the However, the measured suppression seems to be stronger than our expectation, 
indicating that there may be room for other effects. The NA38 group has published @, |J 
also four points for the relative ^' to \1/ suppression in S — U collisions as a function of the 
transverse energy Et- The suppression is larger for larger values of Et, corresponding to 



more central collisions. On the basis of our model we expect such a behaviour, but we lack 
precise quantitative information for the association of a value of Ex to a definite geometric 
configuration. 

In Fig. 3 we show the ty'/ty suppression as a function of xf for Au — Au collisions 
calculated for RHIC and LHC energies. The result is predicted to be energy independent 
at high energies and the two curves coincide. 

The model of coupled channels presented in this letter "naturally" explains, why at 
high energies of the charmonia $ and should be similarly suppressed as observed in 
proton-nucleus collisions at 800 GeV, and why in nucleus-nucleus collisions one should see 
significant differences in and \I/ suppressions as was reported at 200 GeV/ A. 

While the model is free of adjustable parameters, the precision of the two-channel ap- 
proximation is questionable. In order to evaluate the corrections to SV/* from inclusion 
of higher charmonium excitations we switch from the hadronic basis to the quark one in 
coordinate representation. In this case the nuclear suppression of \1/ production is calculated 



as 11 



npA 



(n\v(b, z ,r T )\$f n y 2 



and the same for Here the evolution operator V at high energy (when the fluctuations in 
r T are frozen by Lorentz time dilation) reads, V(b, z, r T ) = exp[— Cr T /2 J 2 °° dz' p^ip, z')}. The 
averaging (...)a in eq. flTT| ) denotes the integration over the coordinates of the production 
point weighted with nuclear density. 

Eq. (JTip (valid only for q = 0) generalizes the two-channel approach in that may 
have other components in addition to the \I/ and ^' states (see an alternative interpretation 



of enhancement of the ty' production rate on nuclei in |TT|, |12||). We have evaluated the 
suppression eq. (|TT|) for \1> and and the relative suppression for various trial 

functions for (r|$^) like exp(— or 2 ), r\ exp(— (3r 2 ), r T K (\r T ) (K is a modified Bessel 
function), where the constants a, f3, X have been chosen to give the same amplitude ratio 
R for the ^' content of |$f^) relative to the \l/ one. The resulting relative ty'/ty nuclear 
suppression exceeds the prediction of the two-channel model 5*^^ = 1 only by 5 — 10% 



for heavy nuclei, a deviation which is still compatible with the data. Another correction 
of the same order of magnitude is expected arising from the x-component of the initial cc 
state. A more complete analysis of these effects as well as recalculation of the nuclear 
suppression using path-integral methods [llj will be presented in a forthcoming paper. 



Note that the initially produced cc wave packet, which is a combination of \I/ and ty', 
attenuates in the nucleus less than each of two charmonia. This fact is important for 
the nuclear suppression of taken separately. This should be checked with available 
experimental data. 

In the case of nucleus-nucleus collisions, we have assumed that charmonia attenuate 
only due to interaction with the projectile nucleons, having xp = ±1. However, particle 
production (and possibly a dense hadronic gas or a quark-qluon plasma) should cause an 
additional suppression of S^T/^, down from our prediction and may explain the deviation of 
our calculations from the results of the NA38 experiment depicted in Fig. 2. Recently, Satz 
13| has proposed to use the value of SV/* as a probe for dense matter formation. 
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Fig. 1 The relative ^'/^ nuclear suppression in p — W collisions. The full circles and 
the solid curve are the data of the E772 experiment M at 800 GeV and our calculation, 
respectively. The open circle and the dashed curve are the result of the NA38 experiment 
|4j at 200 GeV and our prediction, respectively. 

Fig. 2 The relative ^'/^ nuclear suppression in S — U collisions at 200/Ge^. The 
dashed curves are the relative suppression in p — U and S — p (inverse kinematics) 

collisions at 200 GeV . The solid curve, which is the two dashed lines describes our prediction 
for S — U collisions. The data-point is the result of NA38 experiment [Q. 

Fig. 3 Prediction for the relative ^'/^ suppression in Au — Au collisions for the expected 
energies of the RHIC and LHC accelerators. 



